Generalized Quark Transversity Distribution of the Pion in Chiral Quark Models 



O 
(N 



00 
(N 



Alexander E. Dorokhov/'^ Wojciech Broniowski,^'"' and Enrique Ruiz Arriola^'^ 

^ Joint Institute for Nuclear Research, Bogoliuhov Laboratory of Theoretical Physics, 114980, Duhna, Russia 

^Institute for Theoretical Problems of Microphysics, 

Moscow State University, RU-119899, Moscow, RussicQ 

The H. Niewodniczanski Institute of Nuclear Physics, 

Polish Academy of Sciences, PL-31342 Krakow, Poland 

'''Institute of Physics, Jan Kochanowski University, PL-25406 Kielce, Poland 

^ Departamento de Fisica Atomica, Molecular y Nuclear, 

Universidad de Granada, E-18071 Granada, Spain 

^Instituto Carlos I de Fisica Tedrica y Computacional, 

Universidad de Granada, E-18071 Granada, Spaii^ 

(Dated: July 29, 2011) 

The transversity generalized parton distributions (tGPDs) of the the pion, involving matrix el- 
ements of the tensor bilocal quark current, are analyzed in chiral quark models. We apply the 
nonlocal chiral models involving a momentum-dependent quark mass, as well as the local Nambu- 
Jona-Lasinio with the Pauli-Villars regularization to calculate the pion tGPDs, as well as related 
quantities following from restrained kinematics, evaluation of moments, or taking the Fourier- Bessel 
transforms to the impact-parameter space. The obtained distributions satisfy the formal require- 
ments, such as proper support and polynomiality, following from Lorentz covariance. We carry out 
the leading-order QCD evolution from the low quark-model scale to higher lattice scales, applying 
the method of Kivel and Mankiewicz. We evaluate several lowest-order generalized transversity form 
factors, accessible from the recent lattice QCD calculations. These form factors, after evolution, 
agree properly with the lattice data, in support of the fact that the spontaneously broken chiral 
symmetry is the key element also in the evaluation of the transversity observables. 
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I. INTRODUCTION 

The underlying spin-i partonic structure of hadrons 
became first nianifest in the analysis of the deep inelas- 
tic scattering 1]. Actually, further understanding of the 
partonic spin distributions can be gained by the study 
of the transversity distributions [^. From this viewpoint, 
generalized parton distributions (GPDs) jj-Q (for exten- 
sive reviews see, e.g., p-S] and references therein) en- 
code a detailed information on the parton structure of 
hadrons when analyzed at short distances. In the impact- 
parameter space, the GPD's can be viewed as partonic 
probabilities in the infinite-momentum frame distributed 
along the longitudinal momentum fraction (Bjorken-x) 
and the transverse space directions [^, [l3| . It should be 
noted that both GPD's as well as their partonic interpre- 
tation depend strongly on the renormalization scale and 
it is not obvious a priori what, if any, is the reference 
scale, which might have some universal value and signif- 
icance. From a dynamical point of view, the choice of 
such a scale is crucial, as the high-energy modes are inte- 
grated out in favor of an effective and yet unknown non- 
perturbative low-energy dynamics. The renormalization 
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group deals with the intertwining of scales in principle, 
although in practice it can be explored only at the lowest 
orders of the perturbation theory in the running strong 
coupling constant. In addition, GPD's depend also on 
the factorization scheme corresponding to the physical 
process used to extract the partonic distributions at high 
energies. 

/^From a purely theoretical point of view, the great dif- 
ficulty to determine the GPDs from first principles in 
QCD is related to their genuine Minkowski-space nature, 
suggesting application of the light-cone kinematics and 
non-perturbativel y m otivated approaches, such as the 
transverse lattice [ll| , which so far has produced encour- 
aging but scarce results. More recently, however, the 
lowest Bjorken-x moments of the kinematically intricate 
GPDs, the so-called Generalized Form Factors (GFFs), 
have become directly accessible to Euclidean lattices in 
QCD at sufficiently short-distance resolution scales (see, 
e-g-, CI [II)- This is due to the fact that GFFs for 
space-like momenta can be written as matrix elements of 
local operators which can be directly extracted from the 
asymptotics of the Euclidean correlation functions. As 
a further simplification, the scale dependence of GFFs 
in the space-like region undergoes a triangular-matrix 
multiplicative renormalization, which can be easily im- 
plemented (see, e.g., [IJ])- A well known feature of the 
QCD evolution is the loss of resolution at higher energies, 
a property triggered by the existence of the asymptotic 



ultraviolet fixed point, which enhances similarity at in- 
creasingly high Q^ -values. 

In this paper we analyze the quark transversity gener- 
alized parton distribution of the pion (tGPD), related to 
the matrix elements of the bilocal tensor current opera- 
tor q{x)(jfj,,^q{0) (see Sec. ITTl and Refs. [ll, [la| for precise 
definitions). The transversity distribution, also termed 
the maximal helicity GPD, as it involves aligned parton- 
helicity operators, provides insight into the nontrivial 
spin structure of the hadron. For the spin-0 hadrons, 
tGPDs arise due to a nonzero orbital angular momen- 
tum between the initial and final state, and thus offer 
a unique opportunity to learn about the spin structure 
without the many complications of the hadronic spin de- 
grees of freedom, as is the case of the nucleon. Due to 
their inherent complexity, tGPDs are the least investi- 
gated among the hadronic GPD's. In this regard the 
study of the spin structure of the pion is particularly ap- 
pealing and challenging, although at present it is unclear 
how it can be reliably extracted from the high-energy 
experiments. 

The recent lattice determination of the first two X- 
moments of the pion tGPD, denoted as transversity gen- 
eralized form factors (tGFFs) Tr], provides first impor- 
tant and non-trivial information on this issue. The calcu- 
lation was carried out at a lattice spacing of a ^ 0.1 fm 
and a pion mass m^ ~ 600 MeV. For such a small lattice 
spacing the matching to the perturbative MS scheme be- 
comes feasible and corresponds to the scale /i ~ 2 GeV. 
This lattice calculation has triggered some related stud- 
ies focusing either on perturbative aspects of the high- 
Q"^ dependence of the transversity form factors [iq, or 
non-perturbative issues studied within chiral quark mod- 
els [iillOi. 

In this work we analyze the tGPD and the tGFFs of 
the pion for several chiral quark models, extending the 
results presented previously [19[ and providing further 
details. While this unavoidably makes the paper a bit 
technical, we hope that many of the details provided here 
show how a proper implementation of the chiral sym- 
metry, relativity, and normalization can be achieved in 
a non-perturbative model calculation. This is particu- 
larly interesting for the case of nonlocal models, where 
the mass function depends on the momentum. Although 
such models are expected to feature chiral quark dynam- 
ics more realistically, many complications arise due to 
the time-like kinematics implied by the very definition 
of the GPDs. We recall that we are effectively carry- 
ing out the one-loop calculations, where some variables 
are integrated out and some may be left unintegrated. 
Thus, special attention must be paid to the treatment 
of the integrals, particularly to keep the Poincare invari- 
ance explicitly at any step of the calculation, such that 
all results are mutually consistent. 

Via sum rules, the (generalized) form factors are re- 
lated to the GPDs [i-i, mHll]. Experimentally, the 
GPDs of the pion constitute rather elusive quantities 
which appear in rare exclusive processes, such as the 



deeply virtual Compton scattering (DVCS) or the hard 
electro-production of mesons (HMP). 

Chiral quark models have proved to correctly de- 
scribe numerous features related to the vector GPD of 
pion. The parton distribution functions (PDF) have been 
evaluated in the Nambu-Jona-Lasinio (NJL) model in 
Refs. [27l - l29l |. The extension to diagonal GPDs in the 
impact parameter space was carried out in [30.] . Other 
analyses of the pionic GPDs and PDFs were performed in 
nonlocal chiral quark models |3ll - [39| . in the NJL model 
i32i. i40 h 43| I and in the light-front constituent quark mod- 
els [4J, [43. The parton distribution amplitudes, related 
to the GPD via a low-energy theorem |4g] , were evaluated 
in [47 454 1 . The gravitational form factors were computed 
in [55|. Finally, the pion-photon transition distribution 
amplitudes |56l - [59| were obtained in Refs. [6fl-[6J|. 

Besides the phenomenological motivation, it is useful 
to review shortly what aspects of the present investiga- 
tion suggest the use of chiral quark models within the 
present context (see, e.g., [5J]). Firstly, the pion, treated 
as a composite qq state, becomes a Goldstone boson 
of the spontaneously broken chiral symmetry. This of 
course requires the correct implementation of the chiral 
Ward-Takahashi identities - a rather non-trivial point, 
since this condition is not automatically fulfilled in loop 
calculations. At the quark level, this feature is com- 
patible with the large- A^c scaling relations. Within such 
a scheme the pion loop corrections are l/A^c-suppressed 
but chiral-log enhanced at small pion masses. However, 
the leading- A^c contributions present a much milder pion- 
mass dependence, a favorable situation for the unphysi- 
cally large pion masses used on the lattice [iT] . Moreover, 
relativity for the GPDs is properly implemented through 
the so-called polynomiality conditions, and, more specif- 
ically, by the explicit use of the double distributions 
(DDs). Finally, the scale at which a quark model cal- 
culation is carried out can only be identified after a cor- 
rect separation of the momentum fraction carried by the 
quark degrees of freedom. As mentioned already, the par- 
tonic properties depend on the renormalization scale, and 
according to phenomenology |65l . Kjq as well as indepen- 
dent lattice calculations [67[, the (valence) quarks carry 
about 40% of the total momentum at the scale /i — 2GeV. 
In effective quark models, where the quarks carry 100% 
of the total momentum, the perturbative scale is unex- 
pectedly and rather uncomfortably low. However, the as- 
sumption has been tested to higher orders and confronted 
by comparing to a variety of high-energy data or lattice 
calculations. In the present calculation of the transver- 
sity form factors we find again agreement with the data 
after the QCD evolution scheme is implemented, starting 
from a low quark-model scale. 

GPDs in general, and tGPDs in particular, are sub- 
jected to a set of conditions a priori imposed by sym- 
metries and/or completeness, namely, the chiral sym- 
metry, relativity, positivity, and finiteness of sum rules. 
Within the framework of low energy chiral quark models, 
where there is an inherent cut-off marking the low energy 



regime, these conditions are actually not easy to fulfill on 
purely mathematical grounds. Indeed, one-loop integrals 
are four dimensional, whereas GPDs leave two integra- 
tion variables unintegrated and hence some consistency 
is required. However, once this difficulty is mastered, 
which is the case of our approach, there is a trend to in- 
dependence to details of the model. This independence 
is largely enhanced after the QCD evolution, since differ- 
ences are washed out at increasingly higher energy scales. 
The feature is also observed in the study of transversity, 
as to make differences between various chiral quark mod- 
els rather small. 

We apply the local NJL model with the Pauli-Villars 
regularization, as well as two variants of the nonlocal 
chiral quark models inspired by the nontrivial structure 
of the QCD vacuum [63, |63- These models provide the 
results at the quark-model scale. After the necessary 
(multiplicative) QCD evolution liS"], our model results 
are in a quite remarkable agreement with the lattice data 
for tGFFs. Lower values of the constituent quark mass, 
~ 250 MeV, are preferred. 

The outline of the paper is as follows: In Sec. [Tl] we 
give the general definitions of the pion tGPD and tGFFs. 
Then we derive these quantities in the nonlocal chiral 
quark models from the triangle diagram in Sec. IIIIl By 
using the extremely convenient a-representation, we ob- 
tain the corresponding expressions for the tGFFs in the 
momentum- and impact-parameter spaces, the tGPDs 
for the isosinglet and isovector channels, and also, in spe- 
cial forward and symmetric kinematics, the distribution 
of the transversity size of the pion. The analysis is car- 
ried out explicitly for specific nonlocal models in Sec. lIVI 
For numerical estimates of these quantities we use two 
variants of the chiral nonlocal models and the local NJL 
model. In Sec. |V] we present the QCD evolution of the 
above quantities in general, as well as show its conse- 
quences for the studied models. Numerical results for 
the transversity distribution functions after evolution are 
shown in Sec. I VI I Finally, in Sec. lVIII we draw our main 
conclusions. 



II. BASIC DEFINITIONS OF THE 

TRANSVERSITY FORM FACTORS AND 

GENERALIZED PARTON DISTRIBUTION 

In this section we provide the basic definitions as well 
as the kinematics of the transversity observables analyzed 
in the present work. 

The pion u-quark tGFFs, B^'^^ (i), parametrize the 
matrix element 



where the local tensor quark operator is 
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i — y i — ^ 
with D ^ = d ^ — igA^ being the QCD covariant deriva- 
tive, and d ^ = \{ d^ — 9 '^ j . In Eq. ([T]), p' and p are 

the initial and final pion momenta, while P = 2(p' ~^P^' 
q — p' — p, and t = —q^. The symbol TAS denotes sym- 
metrization (S) in i^, /ii, . . . ,/i„-i, followed by antisym- 
metrization {A) in /i, v, with the additional prescription 
that the traces in all index pairs are subtracted (T). The 
factor l/rriTr is introduced b y co nvention in order to have 
dimensionless form factors [17[. Also, as in [l3|, we use 
the positively charged pion and the up-quark density for 
definiteness. 

The above definition, which projects on twist-2 oper- 
ators, can be implemented in a simple and manifestly 
covariant way (see, e.g., [18|) by a contraction with 
two constant auxiliary four-vectors, a and b, satisfying 
a2 = (ab) = and 6^ ^ o. The tGFFs are then defined 
via 
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= (^7r+ [p') 
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u{0)ia'"'a^,b^(iDa] u (0) 
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-1 [(gp) jbp')] V-^ .... yr,U 

2^ (20 ^Tn^ 



i=0, 
even 



where the skewness parameter is defined aq^ 

^ 2 (aP) ' ^ ' 

£, £ [0,1], and (aq), etc., denote the scalar products 
of four- vectors. In Eq. ([3]), [...] denotes the antisym- 
metrization in a and b. 

The tGFFs defined in ([3]) refer to the u-quarks; those 
for the d-quarks follow from the isospin symmetry and 
read 



B^ 



(t)^(-iri?j;:;(i). 



The definition of the corresponding tGPD is [6,] 

{n+ip') I u{~a)ia^'''a^b,uia) \ Tr+{p)) 
[iap)ibp')] 



t-rr J — 1 
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where we do not display explicitly the gauge link fac- 
tor. The tGFFs can be written as the Mellin moments 
of tGPD of the pion as 
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^ Throughout this work we use the so-called symmetric notation. 



III. CHIRAL QUARK MODELS 

In this section we review the generic one-loop features 
of chiral quark models, where the quark self-energy as 
well as the interaction vertices are assumed to have a 
fairly general momentum dependence to be specified later 
on. We derive general expressions for the tGPD at the 
one-quark-loop level, applicable to both nonlocal and lo- 
cal models. We also display formal properties of tGPD 
in our aproach. 



A. Nonlocal chiral quark models 

In the quark-model calculation in the large-A'c limit 
the matrix element (jH]) is given by the triangle diagram 
shown in Fig. [J|. To calculate this diagram we explore 
the manifestly covariant method based on the effective 
approach to nonperturbative QCD dynamics. All expres- 
sions will be computed in the Euclidean space, appropri- 
ate for the process under consideration and, in general, 
for the treatment of nonperturbative physics. The non- 
perturbative quark propagator, dressed by the interac- 
tion with the QCD vacuum, is assumed to have the form 



S{k) 



k + m (fc^) 
D(fc2) 



(8) 



The main requirement imposed on the quark propagator 
is that at large quark virtualities one recovers the per- 
turbative limit. 



s (k) " 4°° 



p- 



(9) 



It is also assumed that the dynamical quark mass, m{k'^), 
is a function rapidly dropping with the quark virtuality 
fc^. It is normalized at zero as 



(10) 



(11) 



m(0) = M„ D(0) = M2. 

We also need the quark-pion vertejo 

Vl{k,q) = ^^,r'^F{kl,kl), 

Jit 

where k± — k zL q/2. The nonlocal vertex F (/c^jfc^) 
is a symmetric function of its arguments, normalized to 
F (fc^, fc^) = m (k^)- In the present study, the nonlocal 
model calculations are performed in the strict chiral limit, 
which means that m (k^ — >■ oo) = 0. 



We should emphasize at this point that the tensor matrix ele- 
ment ^ can not be induced by tadpole- type of diagrams. This 
is evident, because these diagrams depend only on one external 
vector q from which it is impossible to construct the antisymmet- 
ric combination involving the matrix element ((Sjl. In this aspect, 
the results obtained in |20| can not be correct. 
In this work we use the dominant (in the spontaneous symmetry- 
breaking mechanism) structures for the quark propagator and 
the quark-pion vertex. More general structures are used in the 
Schwinger-Dyson approach |70|| . 




FIG. 1. (Color online) The leading-A''c one-quark-loop trian- 
gle diagram contribution to the leading twist tGPD of the 
pion. 



B. Calculation of the triangle diagram 

Within the described approach the triangle diagram 
for the matrix element ^ yields 



MmiU) 



Nr. 



^^^fl 



d^k 



F {k\,k'i) F {kl.k'i) (12) 



-rr{S'(A:+)75S'(fc_)75S'(fc3)cr^^} ( — ^^ — -.a\ af,b^, 

where /c-i- = fc is the initial momentum of the struck 
quark, fcs = fc+ + g is its final momentum, fc_ = fc+ — p is 
the momentum of the spectator quark (cf. Fig. [T|), and 
the covariant average momentum (fc_|_-|-fc3)/2 corresponds 
to the derivative in the definition (|3]). 
After taking the trace one has 
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N, f d^k F {kl,e_) F {klkt) 
4^7 ^ D {kl) D {e_) D [kD 
X {m {k\) [{k^a) (ksb)] - m (fc^ ) [(k+a) (ksb)] 

-m{kl)[{k+a){k^b)]}('^±±^,a'^ ' 



(13) 



where the antisymmetrization in a and b is implied. Con- 
sidering the crossed channel it is easy to get the relation 



{...}('^i±^,a 



(14) 



(-1)" {■••} 



k+ + fc3 



rf— channel 

-, a 



u— cnannei 



in agreement with ([5|). 

For the further analysis, it is very convenient to trans- 
form the integral in (|13p into the a-representation (see 
[ZlllZ2i)i which is one of the basic methods for the study 
of hard processes in perturbative QCD 73], as well as 
in nonperturbative quark models (SJ]. The technical ad- 
vantage of this method is the explicit maintenance of the 
Lorentz covariance. 

Let us define for any function F of virtuality fc^, de- 
caying at large virtuality as 1/fc^ or faster, its a repre- 



sentation (Laplace transform) 



F k')^ 



dat 



-ak^ 



fia) 



(15) 



where F (fc^) is the image of the original function / (a). 
We will use the short-hand F (k^) ~ / (a) . Let us intro- 
duce the following notation [TJ, UM 



2 1,2 



F{kl,kl)F{klk 
D{k\) D {kt) D {kl) 
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D{kl)D{k^_)D{ki) 



m{k1) -- Gm,Q,o{a,(3n), 
m (kl,) r^ Go,m,o (ck, 13, 7) , 
■m{kl) -- Go,o,m(a,^,7): 



(16) 



where the triple a representation (i.e. in parameters a, 
/3, and 7) is applied (see Fig. [T|). With this notation the 
momentum integral in Eq. (J13p is transformed into the 
a-representation expression for the matrix element, 
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N, 
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A3 " V A 

[aGm, 0,0 (a, /3, 7) +/5G'o,m,o (a, /3, 7) +7^*0,0, m (a, /3, 7)] : 
where A = a + /3 4^ 7 and 



d(a/37). 



da / d/3 / d7. 
Jo Jo 



(18) 



The only dependence on ^ in Eq. (J17p appears in the 
polynomial factor in the second line. It is clear that in the 
expansion of this polynomial in powers of ^ only the even 
powers survive, in accordance with Eq. ([3|), since for the 
odd powers of ^ the integrand is antisymmetric in a and 
7. Thus the polynomiality property of Eq. ([7]), namely 
that the X""^ moment of E^ {X, ^, t) is a polynomial in 
^ of the order not higher than n, is immediately evident 
within our approach. 



C. Trans versity pion form factors in momentum- 
and impact-parameter spaces 



^From representation (fT7)) . by using the definition of 
the tGFFs ©, one get^ 



B^mit) 



N^ (n-1)! /d(a/37)^_^, ^^^^ 



ATT^f^i\{n^l-iy.J A 



vn+2 



[2aG™,o,o («, /3, 7) + /5Go,,„,o (a, /?, 7)] /3""'"' ' ^ " 



In the following we will explore the strict chiral limit of m,r = 0. 



where i = 0,2,... < n — 1, and the symmetry proper- 
ties under the interchange of a and 7 has been used. 
The transverse (impact parameter) space representation 
is obtained, by definition, after a 2D Fourier-Bessel trans- 
formation. 



F{bl) 



{2nf 



-'(''^1^)F{t^-ql). (20) 



We then get for even i the expression 
Tm \ ±) IGTT^f^ i\{n-l- z)! J a7A"+i 



-i-i / 7 - a 



(21) 



[2aG„,o,o (a, P, 7) + /3Go,,„,o {a, /3, 7)] /?'' 



D. Pion transversity Generalized Parton 
Distribution 



Through the use of the definition of the tGPD in 
Eq. ([7]) we arrive at the formula 



E^iX,^,t) 



Nr. 



d{aM)^^s^t^ 



(22) 



47r2/2 J A3 

[aGmfl.o (a, P, 7) + /3Go.m,o (a, /3, 7) + 7^0, o,m (a, /?, 7)] 



-l<X.^±i2^<L 

A 

Let us integrate over the /3 parameter, corresponding to 
the quark spectator. From the 5 function we resolve /3 as 



fi^ 



{X + i)a + {X-0^ 
I- X 



(23) 



and apply the positivity conditions for a, j3, and 7. At 
fixed ^ e [0, 1] and X e [—1, 1] one has 3 distinct regions: 

L ^<X<1, where X-hC>0,X-$>0, 
II. - ^ < X < C, where X-K^>0,X-^<0, 
IIL -1<X<-^, where X + ^<0,X-^<0. 

In region I /3 is positive without any limitations. In re- 
gion III all coefficients in Eq. (P5|) are negative, hence the 
support of the integrand has zero measure and the inte- 
gral in Eq. (1^^ equals zero. In the central region II the 
coefficient of a in Eq. (|23p is positive and the coefficient 
of 7 is negative, thus one has the limitation a > 7 ^~y . 
Finally, the total result may be combined as 



E^iX,U)^OiX + 
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47r2/^ Jo ymax{0,7fT#} 

aGm,o,o (a, /3, 7) + l3Go^m,o («, /3, 7) + 7^0,0,™ («, P, 7) 



A2(l-X) 



(24) 



where (x) is the step function, /? is given by Eq. (l23l) . 
and A = [a + 7 + C (a - 7)] /(I - X). 

The isovector and isosinglet tGPDs of the pion are 
obtained as the symmetric and antisymmetric combina- 
tions, 

= E^{X,^,Q^)~E^{-X,^,Q^). 

(25) 

The support of _Ej' ~' is— 1<X<1. The significance 
of the isospin combinations comes from the fact that they 
evolve autonomously with the renormalization scale, see 
Sec. El 



E. Special kinematics: ^ = and £, = X cases 

Some special kinematics is evident. For the case ^ = 
(tPDF) we have 

N r°° 

E^{X,^ = 0,t) = eiX)^-^J d(a7)e-^*x 
2aGm,o,o (a, /3, 7) + f^Go^rn.o (a, P, 7) 



A2(l-X) 



(26) 



Eq. ([22]) and obtain 

pr POO pI 

E^ (X,e,i) = -^ / dL / dxidx2dx3e~='^^^'*x 
47r /^ Jo Jo 

S {1 — xi — X2 — X3) 6 {x — X2 ~ {x3 — a;i) ^) x 

[xiGm,O.Q {xiL, X2L, X3L) + X2Go^rn.O {^iL, X2L, X3L) 
+X3Gofl,m{xiL,X2L,X3L)]. (28) 

To recover the DD representation we further make the 
replacement X2 = b,X3 ~ xi = a and arrive at 

E^ {X, (,t)= db daS {X -b- a^) /| (a, b, t) , 

Jo J-l+b 

(29) 
with the DD identified as 
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^-XiX3t^ 



(30) 



[a:;iGm,o,o (xiL, bL, X3L) + &Go,m,o [xiL, bL, X3L) 

+X3Go,0,m. [xiL, bL, X3L)] . 

Here xi — ^ (1 — 6 — a) and ^3 = i (1 — fe + a). In the 
above expressions the parameter b is non-negative. The 
6 < part of the DD comes from the crossed diagram. 

Sometimes it is also convenient to separate the so- 
called D-term, defined as 



D{b,t) 



1-6 



daf^{a,b,t). 



(31) 



1+& 



where (3 — [a + ^) j^^ ^^^ A = (a -I- 7) jz^- Note that 
in general the first term in the numerator dominates in 
the small X region, while the second one is more impor- 
tant in the region of large X. 

For the border case, ^ = X, we find 



E^{x,^^x,t)^e{x) 



Nr. 



d{aj)e ^*x 



aGm,o,o {a, /3, 7) + PGo,Tn,o (a, /3, 7) -I- 7Go,o,m (a, /?, 7) 



A2(l-X) 



with 13 = 2aj^ and A = [a + J + X {a - j)]j^ 



(27) 



F. Double Distribution 

Some symmetry properties of the GPDs are more 
transparent when they are constructed from the double 
distributions (DDs) [3, [^, UM ■ Actually, the relativistic 
invariance exhibited by the polynomiality conditions is 
manifestly built-in in this approach (see, e.g., Ref. [Z3|)- 
To pass to double distributions, we first make the sub- 
stitution (see, e.g., [761 ) ct = xiL,f3 = X2L,j = X3L in 



G. The b± space and tiie transverse pion size 

Let us now consider tGPD in the transverse coordinate 
space, b±. By using the 2D Fourier-Bessel transform of 
Eq. ((^ one easily gets 

Af i-oo i-co ^2 

E^{X,tbl)=e{X + Oj^ d^ dae-^jy. 
aG„i,o,o (a, P, l) + /3Go,m,o (a, P, l) + 7Go,o,m (a, 13, 7) 



(32) 



Aa7 (1 - X) 



where the value of the parameter (3 is given by Eq. (|23)) 
and A = [a -f 7 + ^ (a - 7)] j^. 

In the zero longitudinal momentum transfer limit, 
^ — >■ 0, one obtains the so-called 3D transverse parton 
distribution 



/?(x,6x) = ii;f (x,e^o,&5 



(33) 



Following [78] one can also introduce the normalized 
quark probability density in the transverse plane. 



PT{X,b^) 



f^{X,b^) 

ffix) ' 



(34) 



where 

f^{X) = E^{X,^ = 0,t = 0), (35) 

as defined in ([25)1 . The partons with the longitudinal 
momentum fraction X occupy within the hadron a disc 
of the average transverse radius squared given by 



bl{X) = j d%^bl}^{X,b, 



(36) 



In chiral quark models the triangle diagram yields 

Ck7 



bl{X)^J^{l-Xf jd{a^) 



(37) 



where /3 = (a + 7) j^y . The C-odd transverse size of the 
hadron, determined by the slope of the tGFF at low mo- 
mentum transfer, can be obtained by integrating b\{X) 
over the momentum fraction, 



dX b\ {X) . 



(38) 



According to Gribov i7i| , one can interpret the normal- 
ized quark density (|M)) as an evolution of the probability 
density for a stochastic motion of a particle in the trans- 
verse plane. The role of the evolution time is played by 
the rapidity variable, 77 = \x\{l/X). For the stochastic 
process one can introduce the mean squared distance of 
the particle as follows j78l |: 



dl{X) = j d%i_blp{X,bA_) = 






(39) 



By using a model with short-range interactions, Gribov 
predicted that ^72] 

dl (77) = Dt^, (40) 

where D is a constant, while in [7g| the result is 



dliv) 



1 



.(1-'^)'; 



(41) 



(47r/.)^ 

Here lo « 0.5 is the slope of the forward quark distribu- 
tion at small X, i.e., q{X) ^ 1/X". Note that Eq. gl]) 
is O(iV-i), since /^ = 0(\/7Vc). Actually, the "chiral in- 
flation" discussed in Ref. i78|] is a pion-loop effect, which 
is 1/iVc-suppressed, but at the same time it is chirally 
enhanced as log(m^) for m^ — )• 0, compared to the lead- 
ing one-quark-loop contribution. In the real world with 
Nc = 3 and m^ = 140 MeV the relative chiral contribu- 
tions to the rms radius of the pion are about 20% 80] B 
Of course, the additional inclusion of pion-loops in our 
model would automatically reproduce this universal in- 
flating phenomenon. 



^ Actually, from the relation for the rms radius of the pion found 
in ChPT [13, (r^) = (T^ - l)/{16n^f^), one has the total low 
energy constant l^ = 13.9 ± 1.3, most of which is saturated by 
the p-meson exchange, l^ ~ 17, at the leading order in Nc- Thus, 
the subleading (1/A^c-suppressed) contribution is estimated to be 



AZ5 = k 



■i',' 



• log(m^/m^) 



IV. MODEL RESULTS 

Having derived the general formulas for tGPDs in chi- 
ral quark models from the triangle diagram of Fig. [1] 
we now pass to presenting explicit numerical calcula- 
tions. We start with the nonlocal models. In the present 
work we consider two variants of the quark-pion vertex 
ofEq. dm). 



Fi {kl, kl) = ^m{kl)m{kl), 



(42) 



FnTY{kl,e_)^-[m{kl)+m{e_)], (43) 

where m (fc^) is the momentum-dependent dynamical 
quark mass. The form (|42l) is motivated by the in- 
stanton picture of the QCD vacuum [68] and is labeled 
"instanton", while the form (|43l) . the Holdom-Terning- 
Verbeek (HTV) vertex, comes from the nonlocal chiral 
quark model of Ref. \6§\ . Some relevant differences be- 
tween both prescriptions regarding the proper implemen- 
tation of chiral symmetry are discussed in Ref. [81[ . 
We consider the dynamical quark mass of the form 



m 



{e) = M,f {e) , 



(44) 



and for simplicity take the profile function /(fc^) as a 
Gaussian, 



f{k^ 



-Ak^ 



(45) 



(note that A has the interpretation of the squared inverse 
momentum cut-off). The model contains two parame- 
ters: the dynamical quark mass at zero momentum, Mq, 
and the nonlocality scale, A. For our numerical estimates 
we take one parameter fixed at a physically reasonable 
value, Mq ~ 240 MeV, and then fix A via the pion d ecay 
constant evaluated in the chiral limit, f^ = 84 MeV !8(|. 
The expression for /^ in the instanton model is given by 
the Diakonov-Petrov formula 68] , 






Nc_ 
47r2 



du'i 



m (u) 



[m{u) — um'{u) + u ra' (u)) 



-,1/2 



(46) 



while in the HTV model one has the Pagels-Stokar for- 
mula [69, 82] 



r 



HTV 



47r2 



du'i 



m (u) 



m (u) um' (u) 



The described parameter-fitting procedure yields 



Ai = 0.7 GeV" 



Ahtv = 0.375 GeV" 



1/2 



(47) 



(48) 



For the instanton model, the integrand in Eq. (|17|) and 
the subsequent formulas can be expressed as follows: 



aGm^o.Q (a, 13, 7) + /3Go,m,o (a, /3, 7) + 7^*0, o,m (a, /3, 7) 

/3dy24</2+7dy24<^ (49) 



i> adl/'dldl^/^ 



while for the HTV model one has 
aGm,o,o (a, /3, 7) + /3Go,m,o (a, l3, 7) + 7^0,0,™ (a, ^5, 7) 



(50) 



+ 


(a <- 


^7) 




+ 


/3 (444 + ^0 


44 + d^id'^^ + did^d^) . 


[ei 


e we 


have intro 


duced the short-hand notation 



D{k^ 



(51) 



For the assumed Gaussian form factor (|15|) the d^ func- 
tion at large a ^ A has the following behavior 



— i— M;e-^("-2"^)e (a - 2nA) , 



with 



R{X) = l-4Am2(A), 
where A is the root of the equation 
X + m^ (A) = 0. 



(52) 



(53) 



(54) 



The functions (|52p can also be used as approximants for 
the analytic calculations of the quark distributions in the 
pion. In the momentum representation this simplification 
means that in the denominators of the integrands we ne- 
glect the momentum dependence of the dynamical quark 
mass, as would be the case of the local quark models. 



A. The numerical results for nonlocal models 

In this subsection we present the results for the nonlo- 
cal models. These results are obtained from the formulas 
presented above with the help of numerical integration. 

We start by exploring the i-dependence. In Fig. [5] we 
present the pion w-quark tGFFs in the HTV model and 
in the instanton model. First of all, the increase of the 
indices n or i causes a decrease of the form factor normal- 
ization. We also note a faster fall-off with t of the tGFFs 
for the case of the instanton model compared to the HTV 
case. We note that the tGFFs undergo the QCD evolu- 
tion, which will be discussed in detail in Sec.|Vl The S^q" 
form factors, however, evolve multiplicatively, hence we 
can read off their i-dependence from Fig. [2] 

At large t the B^q" form factor in the HTV model has 
the asymptotic behavior ^ \nt/t. This follows from the 
asymptotic formula 



B^^„ {t > A-i 

m? (u) 



HTV 1 Nc 



t 167r2/2 



D{u) 



+2 



du- 



, m(u + v) 

dv^rA ^X 



D{u) Jo D{u + v) 

m (u) m{u + v) 



(55) 



n=1,i=0 
n=2,i=0 
n=3,i=0 
3,i=2 




-t [GeV" 



FIG. 2. (Color online) The tGFFs B^f (t) in the HTV model 
(solid line) and in the instanton model (dashed line) for sev- 
eral lowest values of n and i. The sequence in the legend 
corresponds to the sequence of the curves, from top to bot- 
tom. 



n=1,i=0 
n=2,i=0 




0,0 0,5 



b\ [GeV"'] 



FIG. 3. (Color online) The tGFFs S;[f (fe|,) in the impact 
parameter space in the HTV model (solid line) and in the 
instanton model (dashed line). The sequence in the legend 
corresponds to the sequence of the curves, from top to bottom. 



For the instanton model the fall-off is exponential, since 



b: 



TIO 



t 



it » A-1) ^ 



l-2i 



Nc V^M, 



vaTa^ 



4^2/2 ^(A) 
t 



(56) 



Si(. 



AVXt 



In Fig.|3]we display the tGFFs in the impact-parameter 
space. The information is the same as in Fig. [21 as the 
two figures are simply linked with a Fourier-Bessel trans- 
form. Nevertheless, the different large-t behavior of the 
instanton and HTV models is very vividly seen in the 
small-^T behavior in Fig. [31 

Next, we explore the X dependence in the simplest 



-HVT 



-HVT 



2,0 r 




X 




FIG. 4. (Color online) The tPDF in the HTV model (solid 
line) and in the instanton model (dashed line). 



case of t = and ^ = (tPDF). In Fig. |4] we present 
the results of calculations of the tPDF in the nonlocal 
models (|35l) . We notice a more-less triangular shape for 
both models, with a depletion near X — 0. 

The end-point behavior of these functions can be in- 
ferred from Eq. (|26p by using the approximants (j52p . The 
X — > 1 behavior is governed by the properties of the ac- 
tive dynamical quark, while the X — > behavior is re- 
lated to the spectator quark. For the instanton model the 
endpoint behavior is exponentially suppressed, namely 



fi {X ^ I) ^ {1 - Xf e,,p 
2AA 



2AA 



1-X 



fi (X^O)^ cxp 



X 



(57) 



while for the HVT model one has a power-like behavior 

-•^^^ (X ^ 1) ^ (1 - X) , 

f^^'^ (X ^ 0) - const. (58) 



/^ 



We remark here that the end-point behavior in 
Eqs. (I57I58P is sensitive to the radiative corrections, 
hence it evolves with the scale. 

A similar behavior is obtained for the transverse size 
distribution at t = 0, shown in Fig. [51 namely 



bljiX -> 1)- (1-A:)*exp 



2AA 



1- a: 



(59) 



bli {X ^ 0) 

blnvT {X ^ 1) 
blnvT (X ^ 0) 



^exp 

{1-X 
const. 



2AA 



X 



Next, we present our results for the distribution func- 
tion of the mean square distance. In Fig. |6]we show (i]_ 
as a function of X, while in Fig. [7] we present the same 
quantity as a function of the rapidity variable 77. We also 
compare our results to the calculations of Refs. [731 (C-) 



FIG. 5. (Color online) The distribution function of the trans- 
verse size in the HTV model (solid line) and in the instanton 
model (dashed line). 



-HVT 




X 



FIG. 6. (Color online) The distribution function of the mean 
square distance in the HTV model (solid line) and in the 
instanton model (dashed line), plotted as a function of X. 



and [73 (PPV). In the region of large r], corresponding 
to low X, various model predictions are different. 

Finally, we explore the dependence on ^ and X of the 
pion tGPDs at t = 0. The results are given in Figs. [SI 
and [HI We note the symmetry properties following from 
the definition (|^5|) . We can also see that the curves bend 
near X = ^. 

To summarize the study of this subsection we state 
that the results, apart for mathematically different end- 
point behavior, are qualitatively similar in the two ex- 
plored variants of the nonlocal chiral quark models. 



B. Nambu— Jona-Lasinio model 

We term the usual Nambu- Jona-Lasinio model with 
point-like quark-quark interactions the local NJL model. 
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FIG. 7. (Color online) The distribution function of the mean 
square distance as function of rapidity 77 in the HTV model 
(solid line), in the instanton model (dashed line), in the Gri- 
bov approach [7!l| (G) (dot-dashed line), and in the PPV 
model ^Tgj (dotted line). 
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FIG. 8. (Color online) The pion tCPD for isovector case in the 
HTV model (solid lines) and in the instanton model (dashed 
lines) for several values of ^. 



All formulas for the local model follow from the nonlocal 
expressions given above, with the constant quark mass, 
which formally corresponds to taking the limit A — > 0. 
In addition, a regularization prescription, necessary to 
make the divergent integrals finite, is implemented, as 
discussed below. 

The one-quark-loop action of the NJL model is 



r 



NJL 



iiVcTr log {i^ - MU^ 



rcg 



(60) 



where M is the constituent quark mass generated via the 
spontaneous breaking of the chiral symmetry. 



U — exp(i75(/> • t) 



(61) 



with (p denoting the pion field, while m is the current 
quark mass. We apply the NJL with the Pauli-Villars 



FIG. 9. (Color online) The pion tGPD for isoscalar case in the 
HTV model (solid lines) and in the instanton model (dashed 
lines) for several values of ^. 



regularization in the twice-subtracted version of Refs. 154, 

I — t^ — T 

I83 i. 84] . Variants of chiral quark models differ in the way 

of performing the necessary regularization of the quark 

loop diagrams, which may to some extent influence the 

physical results. 

Here we use the prescription where M^ in the loop 

integral is replaced with the combination M'^ + A'^, where 

in the present context A is the cut-off parameter, and 

then the regularized observable is evaluated according to 

the formula 



Cg = 0(0) - ©(A^) + A^d C'(A2)/dA2 



(62) 



Die-multiplying factor gj = AP / f^ is not regular- 



The 

ized [54l.l83ll84j. 

In the local model it is relatively simple to go beyond 
the chiral limit, hence we do not restrict ourselves to the 
case 771^ = 0. Since the lattice data used in this work are 
actually for m^ = 600 MeV, hence not at all close to the 
chiral limit of ttt, — 0, we need to deal with a situation 
of moderately large pion masses. The prescription to fix 
the model parameters is as follows: the three constants 
A, M , and m are traded for the constituent quark mass, 
M, the pion decay constant /jr, and m^. We assume 
that A depends on M only, and not on m. Constraining 
/tt = 93 MeV (the physical value) and using the given 
value of 771^ leaves us with one free parameter only, M, 
which is taken in the 250 — 300 MeV ball park. 

We recall that the optimum value of M used in chi- 
ral quark models depends of particular observable used 
for the fitting procedure. The application to the p me- 
son suggests M above 777p/2 ~ 400 MeV, while the soli- 
ton models for the nucleon prefer M ~ 300 — 350 MeV 
[85j . However, significantly lower values follow from other 
studies in the pion sector. The charge radius of the pion 
in the NJL model with the Pauli-Villars regulator favors 
M ^ 280 MeV 54], however, the pion-loop corrections to 
this observable are important. The analysis of the radii 
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of the pion charge and transitio n for m factors from quark 
triangle diagrams yields M = ^/2jiiTf^ - 240 MeV 3^]. 
Another restriction on the value of M follows from the 
Adler function and the corresponding vacuum polariza- 
tion contribution to the gyromagnetic factor g —2 oi the 
muon. The loo p ap proach (without and with radiative 
corrections) [13, ill yields M = 180 - 200 MeV, the an- 
alytic perturbation model [89| gives 240 MeV, while the 
nonlocal chiral quark model [90] suggests 250 MeV. Our 
chosen value of ~ 250 MeV falls into this ball park. 

In the NJL model the formulas for the lowest two 
transversity form factors are very simple, 



^Tio(^) 

TOtt 



K = 



1 /•! — a 

da I dl3K, 

"'0 
1 /•! — a 

da / d/S aK, 

^0 

NcglM 



27r2 {AP + ml{a - \)a + t^(a + /? - 1)) 



rcg 



(63) 



with ^TT = ^^ / f-n- The variables a and /? are the Feynman 
parameters. 

The result for the tGPD are particularly simple at 
t — Q and in the chiral limit, namely trapezoidal for the 
symmetric (/ = 1) combination. 



Sj'^(X,e,i = 0;Aio)/A^ = 



1, o<x<e ,... 



and triangular for the antisymmetric (/ — 0) combina- 
tion, 

E^^^^iX, e, t = 0; ^io)|N - I ^; ^ |J I 5 . (65) 

Here N denotes a normalization constant following from 
the model. 

Other results of the local NJL model, the correspond- 
ing plots, and comparisons to the predictions of the non- 
local models will be presented in the following parts, to- 
gether with the discussion of the QCD evolution. 



V. QCD EVOLUTION 

We now come to a very important aspect of our anal- 
ysis. Before comparing the results to the lattice data 
we need to carry out the QCD evolution, as the tGPD 
and tGFFs evolve with the scale. The need for the evo- 
lution has been discussed in detail in [43]. In essence, 
our approach consists of 1) evaluation of the appropriate 
soft matrix element in the given model at the low quark- 
model scale, where the matrix element is matched to the 
QCD result, and 2) subsequent evolution to higher scales 
with appropriate perturbative QCD equations. 

For instance, the lattice data correspond typically to 
the scale of about Q = 2 GeV, as follows from the used 



value of the lattice spacing, while the quark model cal- 
culation corresponds to a much lower scale. 



Mo 



A 



QCD- 



(66) 



A detailed discussion of the evolution issue and ways to 
set the quark model scale is presented in Ref. [ij, [91| , 
where the scale 



Aio = 313 MeV 



(67) 



is advocated. We stress that the inclusion of evolution is 
crucial for obtaining the results at experimental or lat- 
tice scales. A non-trivial test is to check that the proce- 
dure reproduces consistently other observables at a given 
scale, ^ (see e.g. Ref. [ij, [91[ for a detailed comparison). 



A. Evolution of tGPD 

The leading-order DGLAP-ERBL evolution for tGPD 
is given, e.g., in [Q]. To carry out this evolution in practi- 
cal terms, we use the method given in |92h95| . where the 
basic objects are the moments in the Gegenbauer poly- 
nomials of index n 



gnifi)^ / dA£;J'^(A,e,i;A*)G3/2(x/0. 



(68) 



The DGLAP region, X > ^, is outside of the orthog- 
onality range for the polynomials Gn {X/S,)- The LO 
DGLAP-ERBL evolution amounts to the multiplication 



g„{lj) = L„5„(/^o), 



J-'n — 



«(m) 

a(/io) 



7X/(2/3o) 



(69) 



(70) 



The anomalous dimensions in the transversity (tensor) 
channel are given by 



7n = y ffn 



8, 



(71) 



where Hn = X]fe=i 1/^- I'^ particular, one has for the 
two lowest form factors 7^ = f and 7^ = 8. We use 
/3o — -tNc — 3 A^/ and the running coupling constant 



a(M)=47r/[/3olog(A.VA?,cD)] 



(72) 



with Aqcd = 226 MeV for N^ 



Nf 



3. The inver- 



sion of the evolved moments back into the evolved GPD, 
applied in our calculation, is explained in |92h95I|. 

We also recall that in the transversity channel the 
quark distributions evolve autonomously, i.e. do not mix 
with the gluon distributions, which is the case of the vec- 
tor and axial channels. That way no gluon tCPDs are 
generated by the QCD evolution, as by construction they 
vanish in chiral quark models at the quark-model scale. 
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B. Evolution of transversity form factors 

The LO DGLAP-ERBL evolution of tGFFs, defined 
as moments of the GPDs, has been spelled out explicitly 
in |14| . The triangular structure which appears from the 
considerations on the evolution of the tGPDs is, for odd 

n = 2k + l, 

k k 

B2k+i,2i = kT{2k) J2 (4"^ + 3)L2™+i ^ (73) 

rn=0 j=k—l 

220-.) (_i)».-.-r {j +m+l) B0,+,,2U-k+i) 
r{2j + l)r(m -j + l)T{k - m + l)T (fc + m + |) ' 

and, for even n = 2fc + 2, 

fc fc 

B2k+2.2l - r(2fc + 2) ^ (4m + 5)L2,n+2 J2 C^^) 



In the chiral limit and at i = 



?n— 



j=k-l 



22,-2fc-l(_^)n^-,r (j + ^ + I) i3,V,) ,(^._,^,) 

r{2j + 2)r(m -j + l)T{k - m + 1)T (fc + m + |) ' 

where fc = 0, 1, 2, . . . and I = 0, 1, . . . , fc. We have in- 
troduced a short-hand notation Bni = B^^^lt; p) and 
Bjjj = i?5^^,-(i; /^o)- For the lowest moments we have, 
explicitly, 

Bio = LiBiQ, 

B32 = Ti^i- L'i)B^Q + £3^32, 

B54 = -^i^L^ - 14^3 + 5i5)B?o 

2 
+ (-^3 - -^5)^32 + L^B^^, 



B20 — L2B20, 
3 

-B42 = y (-^2 — Li)B20 + ^4-^42, 

^30 = -^3^30' 

2 
-B52 = 0(^3 - -^5)-B3o + L^B^2i 



-B40 — ^4^40. (75) 

In particular, the two lowest tGFFs available from the 
lattice data, Byj^ and BI^'^q, evolve multiplicatively as 



follows: 



^Tnoi*'^ ^) = ^Tnoi^i A^o) 



which numerically gives 



q(m) 

a(/^o) 



7^/(2/3o) 



B;fo(i;2 GeV) = 0.75i3j'i"o(t; /io), 

Bj;,';(t; 2 GeV) = 0.43Bj2"o(i; /io). (76) 

Note a stronger reduction for Bt20 compared to Btiq as 
the result of the evolution. 



^Tioi^ = 0;/io)/TO^ 



B^2o{t ^ 0; H) 
■Stio(^ = 0;^) 



q(a^) 

a(^o) 



8/27 



(77) 
(78) 



VI. 



NUMERICAL RESULTS AFTER THE QCD 
EVOLUTION 



In this section we present our numerical results after 
the QCD evolution for the tGPD of the pion, its special 
cases '^ = and ^ = 1, corresponding to the tPDF and 
tDA, respectively, as well as discuss the tGFFs. The 
latter are compared to the available lattice data of Ir] . 



A. tGPD 

The results of the calculation of the tGPD of the pion 
at a sample value of ^ — 1/3 and at i = 0, together with 
the LO DGLAP-ERBL evolution, are given in Figs. [TU] 
andim For the non-local case we take the HTV model 
(I43l50p , as the results of the instanton model (I42I49P are 
qualitatively similar. Here we take for simplicity the chi- 
ral limit, ttIt^ = 0. We provide in the figures the symmet- 
ric (S) and asymmetric (A) combinations in the X vari- 
able (PSJ). The solid lines correspond to the calculation 
at the quark-model scale, ^o- In this case we convention- 
ally normalize the plotted functions with a constant N 
in such a way that 



dXE:^-^{X,C,t 



0;txo)/N- 



1 + e 



(79) 



for all displayed models. 

Further, we note the gross qualitative similarity be- 
tween the nonlocal HTV model and the local NJL model. 
The differences are manifest in the end-point behavior. 
Near X = 1 the tGPD in non-local model is suppressed, 
as explained in Sect. IIV Al Also, near A" = the quan- 
tity E^' is depleted compared to the local case, where 
no minimum is present. 

The dashed and dotted curves show the results evolved 
to the scales 2 GeV and 1 TeV, respectively. After the 
evolution the results of the HTV model and the local NJL 
model are qualitatively very similar. 



B. tPDF 

Next, we explore the special case £, = 0, again for t ^ Q 
and TOtt — 0. In this case tGPD corresponds, by defini- 
tion, to tPDF. In Fig. [12] we compare the predictions of 
the three considered models at the quark-model scale, 
fiQ. We note different end-point behavior, both at X = 1 
and at A" = 0, according to the discussion presented 
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FIG. 10. (Color online) The DGLAP-ERBL evolution of the 
symmetric (S, or / = 1) and antisymmetric (A, or / = 0) parts 
of the quark tGPD of the pion in the non-local HTV model 
for m^ = 0, i = 0, ^ = 1/3, and M = 240 MeV. The solid 
line corresponds to the initial condition at the quark model 
scale fio ~ 313 MeV, the dashed line shows the result of the 
evolution to /x = 2 GeV, and the dotted line to jj, — 1 TeV. 



in Sect. IIV Al Near X — 1 the instanton model has a 
stronger suppression in tPDF than the HTV model. The 
local model approaches zero linearly. Again, we note that 
the QCD evolution changes the end-point behavior. 



C. tDA 

Another interesting limiting case is provided with 
^ = 1. In that case 



E^{X,t^O,^^l)^^T{X), 



(80) 



where (j)^{X) is the tensor distribution amplitude of the 
pion, defined as 



{0\d{z)aa/3j5ui-z)\TT+ (q)) = 
•J Jo 



(81) 



where X = 2u — 1 and N'^ is the normalization factor 
yielding J^ du(l)T{u) — 1. 



FIG. 11. 

model. 



(Color online) Same as Fig. [TO] for the local NJL 



The local NJL model predicts a constant (p^{X) at 
the quark-model scale. Again, as seen from Fig. [131 the 
difference between the local and non-local models is seen 
in the end-point behavior, X '^ ±1. In the intermediate 
range of X the tDA (t)J^{X) is close to a constant also for 
the non-local models. 

In Fig. [H we show the LO ERBL evolution of the tDA 
of the pion in the local NJL model. We note a gradual 
approach towards the asymptotic form 



'^T,asym 



(u) — 6m(1 — u). 



(82) 



For the non-local models the effect of the evolution is 
similar. 



D. tGFFs 

In Fig. [15] we show the LO DGLAP-ERBL evolution 
of the tGFFs evaluated in the local NJL model. 

By comparing the two panels of Fig.[T5|we note that for 
the tGFFs is multiplicative, and increasing the scale leads 
a quenching of Bxno the form factor. For the form factors 
Btui with i 7^ the evolution is more complicated, as can 
be inferred from Eq. (1751) . For the non-local models the 
effects of the evolution for tGFFs are similar. 



14 





1.2 


1 1 


^ 


1.0 


\ 


o 
II 


0.8 
0.6 


/ 
- / 


X 


0.4 


" 1 : 



^ 



0.2 



0.0 




0.0 



0.2 



0.4 0.6 

X 



1.0 



FIG. 12. (Color online) Comparison of the tPDF {E^{X,t = 
0,^ = 0)) in the local model (solid line), instanton model 
(dashed line), and the HTV model (dotted line) for tti^ = 0, 
evaluated at the quark-model scale. 




FIG. 14. (Color online) Evolution of the tensor distribution 
amplitude, tDA, in the local NJL model. The subsequent 
curves (from bottom to top &t u = 1/2) correspond to /i = 
/io = 313 MeV (the constant), ^ = 500 MeV, ^ = 2 GeV, 
/i = 1000 GeV, and fi = oo (the asymptotic form 6ii(l — u). 





'. _ 


', ' 


1.0 


- ^ — — '-^_ — _ 


— — ^ ^ ~~ — ^ \^ 


/; 




^ 0.8 






R '■' 


1. 
i; 


',1 - 
'I - 


$ 0.4 


1' 

- 1,' 

1' 


.1 - 
',1 - 
■ 1 : 


0.2 


1. 
i; 


',1 " 
'1 - 


0.0 


-i 


1^ 



models yield the nornialization 



-1.0 



-0.5 



0.0 
X 



0.5 



1.0 



FIG. 13. (Color online) Comparison of the tDA (0t(X)) in 
the local model (solid line), instanton model (dashed line), 
and the HTV model (dotted line) for m^ — 0, evaluated at 
the quark-model scale. 



E. Chiral quark models vs lattice 



The content of this Section has aheady been presented 
by us in a greater detail in 19]. For the completeness of 
the present work we repeat the main results. 



The presently available fuU-QCD lattice results 117[ are 
for Bfo" and ^26" a-nd for —t up to 2.5 GeV^, with mod- 
erately low, but still away from the physical limit, values 
of the pion mass, m,^ ~ 600 MeV. The calculation of [13] 
uses the same Nf = 2 set of the QCDSF/UKQCD en- 
sembles with improved Wilson fermions and the Wilson 
gauge action that were used previously in the analysis of 
the pion charge and gravitational form factors [96J . 

We note that for t — both the local and non-local 



Brr\n{t — 0;/io)/m7r — ST^/^ 



^ Io° ^^ (uIZhu))-' ^^^^^ ~ w™'(w))' (83) 

-TfTT.u j , n \ / A'"^ I r*^^ J um(u) 

X (to^(u) -I- ^um{u)m'{u) + ^u'^m''^{u)) 
-irduj^^^^{m{u) + 2m^iu)m'iu))}, (84) 

where m'{u) = dm{u)/du. In the local limit, where 
m{k'^) — >■ const, one reproduces Eqs. (177178^ . 

The results for i3j^"o(i)7 ^ = 1; 2, are shown in Fig.lTSl 
In our study we have assumed that BTno/triT^ depends 
weakly on m7r, similarly to the local model case [19|. 
Therefore, to compare to the lattice data for Bmo, we 
multiply the results of the calculations obtained in the 
chiral limit with 171,^ — 600 MeV. We have carried out 
the QCD evolution procedure as described in the pre- 
vious Sections, from the quark model scale up to the 
lattice scale of 2 GeV. From Fig. [11] we note that the 
HTV model with the vertex function given by Eq. p3)) 



(solid lines) and with Mg = 300 MeV works best, de- 
scribing accurately the data, while the instanton model, 
Eq. (|^^ (dashed lines), results in form factors falling-off 
too steeply. We have found that lower values of Mq spoil 
the agreement with the lattice data. 

In Fig.[T7]we show the results from the local NJL model 
evolved to the lattice scale oi ^ = 2 GeV, confronted with 
the lattice data scanned from Fig. 1 of [17]. We have 
used m^ = 600 MeV and selected M = 250 MeV, which 
optimizes the comparison. As we see, the agreement is 
remarkable. 

In Ref. [i9| we have also investigated the dependence 
of the values of the form factors at t = on the value of 
TO^ , as studied in [13] • We have also noted in 19] that the 
results presented in Fig. [T7] depend quite sensitively on 
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FIG. 15. (Color online) The transversity form factors Bni{t), 
evaluated in the local NJL model at the quark-model scale 
fj,o (top panel) and evolved to /i = 2 GeV (bottom panel). 
Solid line - -Bl'o(i), dashed line - B^o{t), dotted line - B^oit), 
dash-dotted line - B^2{t)- 



the value of the constituent quark mass, M, with higher 
M yielding lower values of the transversity form factors. 



VII. CONCLUSIONS 

In the present paper we have shown how the spin- 
less pion acquires a non-trivial spin structure within 
the framework of chiral quark models. This has been 
achieved by computing the transversity distributions, 
corresponding to matrix elements of the tensor quark 
density, within chiral quark models, where the pion arises 
as the pseudo-Goldstone boson of the spontaneously bro- 
ken chiral symmetry. Moreover, we have worked at the 
leading order in the 1/Nc expansion, which amounts to 
carrying out one-quark-loop calculations, where the im- 
plementation of the symmetry constraints becomes ab- 
solutely essential. Chiral symmetry is respected by im- 
plementing the pertinent chiral Ward-Takahashi identi- 
ties at the quark level. Moreover, the relativity con- 
straints are fulfilled in terms of the polynomiality con- 
ditions which are manifestly preserved through the use 




t [GeV 



FIG. 16. (Color online) The transversity form factors in the 
HTV model (solid line) and in the instanton model (dashed 
line). The data come from |9a ]. 




-t [GeV2] 

FIG. 17. (Color online) The transversity form factors ob- 
tained in the NJL model (lines) for M = 250 MeV and 
rmr = 600 MeV, evolved to the lattice scale of 2 GeV and 
compared to the lattice data from Fig. 1 of [l3] (points). 



of the double distributions, or, equivalently, by working 
with the a-representations. 

We have provided comprehensive results for the tGPDs 
of the pion, as well as related quantities following from re- 
strained kinematics, evaluation of moments, or taking the 
Fourier-Bessel transforms to the impact-parameter space. 
We have also shown in detail various technical aspects of 
our analysis, including the use of the a-representation in 
the nonlocal models. 

The generated tGPDs are defined at a given low-energy 
quark-model scale, and comparison to data or lattice 
results corresponds to implementing the suitable QCD 
evolution. Actually, while the momentum-transfer or, 
equivalently, the impact-parameter dependence of the 
tGFFs remains scale independent, their absolute normal- 
ization does depend multiplicatively on the renormaliza- 
tion scale. Remarkably, the absolute predictions for the 
multiplicatively evolved -BtmOj for n = 1,2, agree surpris- 
ingly well with the lattice results, supporting many pre- 
vious calculations following the same chiral-quark-model 
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scheme amended with the subsequent QCD evolution. 
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